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O^ ' We study the long-time behaviour of a nonlinear Fokker-Planck equation, which 

models the evolution of rigid polymers in a given flow, after a closure approximation. 
The aim of this work is twofold: first, we propose a microscopic derivation of the clas- 
sical Doi closure, at the level of the kinetic equation ; second, we prove the convergence 
of the solution to the Fokker-Planck equation to periodic solutions in the long-time 
limit. 



^ ! 1 Introduction 

In a previous work |TTj . two of us have studied the long time behaviour of some flows of 
infinitely dilute flexible polymers. It has been proved that for appropriate boundary data, 
and provided the solution is assumed regular, the solution returns to equilibrium in a long 
time, whatever the initial condition. The mathematical ingredient for the proof is the 
appropriate use of Log-Sobolev inequalities and entropy methods d la Desvillettes- Villani, 
applied to the Fokker-Planck equation derived from statistical mechanics and kinetic theory 
considerations and modelling the evolution of the polymer chains. It turns out that flows of 
rigid polymers exhibit equally interesting and, actually, much more varied properties in the 
long time, in particular stationary periodic-in-time motions. Several behaviours have been 
experimentally observed. Numerical simulations confirm the ability of the models employed 
to reproduce such behaviours. We refer for example to [UlSKT^ for previous mathematical 
studies. Those behaviours are traditionally filed in different categories with an appropriate 
terminology. One speaks of flows exhibiting kayaking, tumbling, etc. Mathematically, 
one underlying question that might be considered is to investigate whether the solution 
to the Fokker-Planck equation ruling the evolution of the microstructure (here, typically 
rigid rods) converges in the long term to a periodic-in-time solution. This is in sharp 



contrast to the case of flexible polymers considered in |11| where the long time limit of the 
solution to the Fokker-Planck equation is a steady, that is time-independent density. The 
mathematical ingredients mentioned above (Log-Sobolev inequalities and entropy methods) 
are again useful, but their use is more delicate, as will be seen below. 

The purpose of this article is to consider a simple setting, where the long term behaviour 
of the evolution of the microstructures can be proven to indeed be periodic. 

To start with, we consider a commonly used model, namely the rigid rod model with a 
Maier-Saupe potential. It has been observed numerically [18J that, for this specific model, 
the flow is, in the long time, periodic-in-time. It is therefore an adequate setting to consider 
with the view to proving that, in the long time, the solution to the Fokker-Planck equation 
converges to a periodic-in-time solution. This Fokker-Planck equation formally writes 

— = LWM/ (1) 

where L(^) is a nonlinear nonlocal partial differential operator, essentially parabolic (see 
the weak formulation ([6]) below) and ^ is the probability density function describing the 
state of the microstructure. We are unfortunately unable to prove mathematically the 
expected long term behaviour of the rigid rod model with a Maier-Saupe potential. Note 
that this is indeed a particularly challenging issue to make a period appear in an equation 
of the form ([1| which does not explicitly contain any periodic function to start with! In 
some sense, we would need a Poincare Bendixson type theorem for an inflnite dimensional 
system. This is beyond our reach. 

So we proceed somewhat differently. In Sections 12.21 and |3l we first derive an ap- 
proximation of the model, that gives rise to a closed evolution equation for the so-called 
conformation tensor. This equation agrees with the equation obtained when a classical 
Doi-type closure is performed on the original model. In passing, we motivate in Section [312] 
our particular choice of approximate model. We finally prove, in Section UJ that the solu- 
tion to this equation does behave as expected in the long time: it becomes periodic-in-time. 
Our proof falls in essentially two steps. We first show, in Section [4.11 that the conforma- 
tion tensor x^x"^ calculated from the solution ^ to our approximate model, and which 

satisfies the closed evolution equation ([8]) (this is the whole point of the closure approxi- 
mation), becomes periodic-in-time in the long term. We next use this result in our final 
Section 14.31 to conclude our study on the convergence of the density ^ itself. Because our 
fundamental tool (in the course of Section 14. ip is the Poincare-Bendixson Theorem (fol- 
lowing the work |14|). our main result is unfortunately restricted to the two-dimensional 
setting. Other intermediate results we prove however hold whatever the dimension (and 
they have been proved and stated so). Another technical limitation lies in the fact we 
exploit, for our proof, the specific explicit expression of the time-periodic solution we es- 
tablish the existence of, and which attracts all solutions in the long-time. More generality 
in the technique of proof would be highly desirable but is out of our reach to date. 

In summary, the main contributions of this work are of two different types in nature: 
from a modelling viewpoint, we propose a microscopic derivation of the quadratic Doi 
closure, using a stochastic dynamics with constraints on an average quantity (see Propo- 
sition I3.ip ; from a mathematical viewpoint, we analyze the longtime behaviour of the 
solution to a nonlinear Fokker-Planck equation which converges to a periodic in time func- 
tion (see Proposition 14. 6p . 

Our study can be considered in the vein of several previous studies such as [U |3] . It is the 
authors' wish that the quite limited study performed here will be yet another incentive for 
mathematicians to consider the question of long term convergence to non steady stationary 
states for solutions to kinetic equations of the type ([T]). We reiterate that, in our opinion. 



the issue of proving, in some particular settings and under appropriate conditions, that 
solutions to nonlinear Fokker-Planck equations of the type ([T|) converge in the long time 
to periodic-in-time solutions is an interesting, unsolved mathematical issue. 

2 The Original model and the Doi closure 

2.1 The Maier-Saupe model for rigid polymers 

Using the notation of |14j, we consider, in Stratonovich form, the following stochastic 
adimensionalized model for a rigid polymer with the Maier-Saupe potential: 

dXt = P{Xt){KXt + mE{Xt Xt)Xt) dt + V2P{Xt) o dBt, (2) 

where Bt denotes a d-dimensional Brownian motion, A^ is a dimensionless concentration 
parameter and the matrix k € M'^^'^ is related to the velocity field of the ambient flow, 
which is assumed to be homogeneous (so that transport terms in ([2]) are omitted). The 
purpose of the projector operator 

pm = id-^ (3) 

||v\ II 

is to ensure the preservation of the rigid polymer (nematic crystalline polymer) norm \\Xt\\ 
in time: 

d\\Xtf = 0. 

Here and in the following, we use the tensor product notation: for two vectors u and v in 
M , u (^ V is the R matrix whose (i, j)-entry is UiVj. Notice that P{X) is a symmetric 
matrix such that P{X)P{X) = P{X). 

The right-hand side of Equation ([2]) contains three terms which model different phe- 
nomena. The first term models the reorientation due to the velocity gradient k of the fluid. 
The second, non-linear term contains the force associated to the Maier-Saupe potential, 
which describes the effective interaction of a rigid rod with the other rods. The rightmost 
term is a rotational diffusion term. 

In the following, we will in particular consider the two-dimensional case (d = 2) and a 
simple shear flow, for which 

_ Pe r a + 1 " 
""' 2 a-l 



(4) 



where Pe (the Peclet number) and a (a molecular shape parameter) are two adimensional 
parameters, see |14j . But for the moment being, we consider ([2]) in full generality. 
Equation ([2|) equivalently writes, using Ito's integration rule: 

dXt = P{Xt){KXt + 4NE{Xt Xt)Xt) dt + V2P{Xt)dBt - {d - 1)^^ dt, (5) 

\\^t\\ 

where d is the dimension of the ambient space. Obtaining ([5]) from ([2]) is straightforward 
using (with implied summation over repeated indices) 

V2PiAXt) o dBl = V2Pi,j{Xt)dB{ + ^{V2Pj^k{Xt)dj) • V2Pi^k{Xt)dt, 



and 



Pj,k{^)9jPi,k{x) = {6j,k - XjXk/lxl"^) dj{5i^k - XiXk/\x\'^) 
= {l-d)xi/\x\\ 



where 6ij is the Kronecker symbol. 

The Fokker-Planck equation associated to ([5|) (or equivalently to ([2])) writes, in weak 
form: for any smooth test function (/9 : M"^ ^ M, 



with again imphed summation over repeated indices and where for any time t >0, fit{dx) 
is the law of Xt (with support on a sphere) . 

2.2 The Doi closure to derive a closed equation for E(Xt (X) Xt) 

As announced in the introduction, we now recall how to derive from ([5]) a closed evolution 
equation on the conformation tensor 

Mit)=K{Xt(E)Xt), 

using the standard quadratic Doi closure approximation [7] . It is the aim of Section 13.21 
below to provide a microscopic justification of this closure. 

Using elementary Ito diff^erential calculus (and the properties P = P^ and P"^ = P), 
we compute: 

d{Xt ® Xt) = dXt ^Xt + Xt^ dXt + 2PP'^iXt)dt 

= {P{Xt)KXt ®Xt + Xt® XtK^P{Xt)) dt 

+ m {P{Xt)E{Xt Xt)Xt ^Xt + Xt^ XtE{Xt <S) Xt)P{Xt)) dt 

- 2{d - l) ^'®^' dt + 2P{Xt) dt 

\\-^t\\ 

+ V2P{Xt)dBt 0Xt + V2Xt dBtP{Xt) 
= {P{Xt)KXt ®Xt + Xt® XtJP{Xt)) dt 

+ m (p(Xj)E(Xi ® Xt)Xt ®Xt + Xt® Xt^Xt ® Xt)P{Xt)) dt 

+ 2fld-d^^V^ 
+ loc.mart. 

where loc.mart. denotes a local martingale that we do not need to make precise for the 
rest of our argument. Taking the trace of the previous equation and using that ii{AB) = 
iY{BA) and {X ® X) P{X) = P{X) {X ® X) = 0, we check the preservation of the norm 

d\\Xtf = ^ 

as was announced earlier. We henceforth set \\Xt\\ = L. 
Taking now the expectation, we obtain: 



at L^ 



kM + Mk^ - -;r7TE [K-.Xt® Xt Xt (g) Xt) 



+ m Um^ - -^E (M ■.Xt(^XtXt(^ Xt) ] (7) 

+ 2Id-||M, 



where we here introduced the Frobenius inner product: for two M.'^^'^ matrices A and B, 
A:B = tT{AB^) = Zt,=iA,BiJ. 

At this stage, we use the so-called quadratic Doi closure |7j that consists in performing 
the following approximation: for any deterministic matrix K, 

E {KXt ®XtXt®Xt)=K: E(Xi ® Xt) E(Xt ® Xt). 

The following closed nonlinear first order differential equation that rules the evolution of M 
in time is thereby obtained (using the fact that tr(M) = L?'): 

dM rr 2 

= kM + Mk^ — -K, : M M 

dt tr(M) 

+ 4n(2M^ ^M:MM] (8) 

V tr(M) ; ^ ' 

+ 2Id-^^M. 
tr(M) 

Note that, at this stage, the above equation is formal since we do not know that tr(M) 
does not vanish. It is a consequence of the following proposition that this is not the case. 

Proposition 2.1 Assume that we supply equation ([8|) with an initial condition M(0) that 
is a symmetric matrix and that satisfies tr(M(0)) = L^ > 0. Then there exists a unique 
solution M{t) to ([8]). Moreover, this solution remains symmetric for all times and with 
constant trace: tr(M(t)) = L^. 

Proof. We consider a time interval on which equation ([8]) is well posed. Such a time 
interval exists by a standard application of the Cauchy-Lipschitz Theorem. Indeed, the 
right-hand side of ([8|) is a rational function in the coefficients of M. Momentarily, this 
time interval may be bounded, but we will soon see that it is in fact infinite. 
The adjoint matrix M^(t) of M{t) then satisfies the following equation: 

kM^ + M^K^ ^K : M M^ 



dt tr(M) 

+ 4N I 2(M^)2 ^M : M M^] 

V tr(M) ; 



2d rj. 

+ 2 Id ttttM^. 

tr(M) 

It follows that, when M{t) solves (j8|, both M(t) and M (t) are solutions to the first 
order evolution equation 

dB rj. 2 

—- = kB + Bk^ — -k: MB 

dt tr M 



j2 



2 



+ 2B^ — -M : MB 

tr(M) 

+ 2Id-^^S, 

tr(M) 

and that this holds for the same initial condition M(0) since the latter is symmetric. 
Now, the right hand side of this differential equation is a second order polynomial in B, 
with coefficients that are obviously continuous in time (in turn because M{t) solves dB])). 
It follows that the Cauchy-Lipschitz theorem holds for this equation, and thus that the 



solution is unique for a given initial condition. This proves that M{t) = M'^{t) for all 
times in the considered time interval. 

We now take the trace of ([8]) and we easily check that 

|tr(AO = 0, 

using the fact that M is now known to be symmetric. The trace of the solution is thus 
preserved in time and this in particular shows that the solution to ([8]) is defined for any 
time (for any initial condition with non zero trace). We finally notice, as this will be useful 
below that, since tr(M)|t=o = L'^j (El) also writes 

dM r 2 

dt L2 

+ 2M2-4tM:MM (9) 

+ 2Id-^M. 





Remark 1 A natural question is whether Equation ([8]) also preserves positiveness (in the 
sense of symmetric matrices). This property will be a consequence of a rewriting of the 
solution to (HI) as M{t) = E(Xt ® Xt) for Xt solution to a modified stochastic differential 
equation, see Section [21 below. We are unable to prove this preservation otherwise. 

2.3 A rewriting of the equations 

It is enlightening to compare our equation ^ with the equation (10) in the article |14| by 
G. Forest and collaborators. For this purpose we recall the three dimensionless numbers 
used by these authors: the molecular shape parameter a, the Peclet number Pe and the 
dimensionless concentration number N. In this context, the dimension is d = 2, the length 
is L = 1 and k is (|4]) and thus writes: 

K = Fe{n + aD) 



where Q. = ^ 



1 
-1 



and D = 2 



1 

1 



The equation ([7|) then reads (using the fact that Q is skew-symmetric): 
dM 



dt 



: Pe {nU -Mn + a{DM + MD) - 2aE {D:Xt® Xt Xt ® Xt)) 

+ 8N (M^ -E{M -.Xt^XtXtC^ Xt)) 
+ 4(Id/2-M). 



(10) 



To agree with the notation of [13j, we introduce Q = M — Id/2, the traceless part of M. 
Equation (|10|) then rewrites: 



dQ 

dt 



: Pe {QQ -Qn + a{DQ + QD + D)- 2aE {D : Xt® Xt Xt ® Xt)) 

+ m (Q2 + Q/2-¥.{Q:Xt® Xt Xt ® Xt)) 
-4Q. 



(11) 



We observe that the above equation agrees with the Equation (10) in |14| up to multiphca- 
tive constants that do not affect our conclusions. Using now the Doi closure approximation, 
we finally get 



dQ 
dt 



Pe 



nQ-Qn + a{DQ + QD + D)- 2aD : Q ( Q + - Id 
Q - 2iV (^Q + i Id j Q + 2NQ : Q (q + ^ Id 



(12) 



3 Derivation of an evolution equation for Xt that yields our 
closed equation on M 

We have derived in the previous section a closed equation ([8]) on M, using a classical closure 
technique [a la Doi) on the original dynamics ^ on X^. The question now naturally arises 
to know whether it is possible to modify the original stochastic dynamics ([5]) itself so that 
M(t) = K{Xf Xf) calculated from Xj solution to this modified dynamics is solution 
to dSD. 



3.1 Two possible closures on the stochastic differential equation 

To begin with, using the fact that for any vector x G M and matrix A G M , (x^x) Ax 
A : (x x) X, we write ([5]) under the form: 

dXt = {KXt + mE{Xt Xt)Xt) dt 

^ (k : (Xt Xt)Xt + ANE{Xt Xt) : {Xt Xt) Xt) dt 



\Xt 



Xt 



+ V2P{Xt)dBt -{d- 1)^^ dt. 

\\^t\\ 

We next modify this equation as follows: 

dXt = {nXt + 4iVE(Xt ® Xt)Xt) dt 

- En^jpy ('^ • ^(^* ® ^*) ^* + 4A^E(^t ® Xt) : E(Xi ® Xt) Xt) dt ^^^^ 

+ V2RtdBt-X—^^-^dt, 

where the pair (Rt, A) is yet to be determined so that M(t) = E(Xt Xt) calculated from 
Xt solution to (|13p is indeed solution to ([8]). Here, Rt G M is an adapted stochastic 



process, and A G M is a deterministic constant. 
As above, an elementary Ito calculation yields 

d{Xt Xt) = dXt ^Xt + Xt^ dXt + 2RtR^dt 

= UXt CdXt + Xt® XtJ + 8iVE(Xt Xt)Xt Cd xi) dt 

{2k : E{Xt ® Xt) Xt®Xt + ^ NE{Xt ® Xt) : E{Xt Xt) Xt Xt) dt 



mxtW) 

Xt ^ Xt rp 



+ loc.mart. 



Taking the expectation, we obtain 



— = (kM + Mk^ + 8NM^) 
dt 

7TTT (2k : MM + 8NM : M M) 

tr(M) ^ ^ 



This equation is then equivalent to equation ([8]) if and only if 



tr(M) ^ ' t ^ tr(M) ' 

Obviously, the simplest possible choice for {Rt,X) is to set 

Rt = Id and X = d. (14) 

It follows that the diffusion term in the associated Fokker-Planck equation is simply a 
Laplacian. Let us write the non-linear Fokker-Planck equation we thus obtain: 

dt -^'""[y ^^+tr(M[V']) J^ 

+ 4iVdiv (^i^-M[^P]x + ^^j|^:|^^ ) V' ) (15) 



M[V'(t, •)]=/ x0xi;{t,x)dx. (16) 



where 



We will see in the next section that this particular choice (|14|) of pair {Rt,X) may be 
motivated by modeling considerations. This turns out to be the choice we advocate. 
An alternate convenient pair (among many possible choices) would be to set 

RtRf = Id -— and X = (d - 1). (17) 

^ * tr(M) ^ ' ^ ^ 

Using the fact that M < tr(M) Id (in the sense of symmetric matrices), the existence of 
such a matrix Rt follows from Cholesky factorization, for example. Then, the associated 
Fokker-Planck equation would write 

dt VV tr(M[V;]) J"^ 

+ AV' - -4^^ : VV + id- l)div ( ,„^,,,, V 
^ tr(M[V']) ^ ^ ^ ytriMl'ip])^ 

where M[ijj] is defined by p6|) . We have not been able to motivate the alternate choice ([Tj 
as convincingly as the choice (|14p . and will show that our preferred choice (J14p enjoys 
several agreeable properties. We will therefore henceforth adopt ([1] 



3.2 A possible justification of our approximation (IT^ 

To derive an appropriate approximation of Equation ([5]), we now follow a different path. 
Since ([5]) is the projection on the manifold defined by the constraint "||Xj|p constant" of 
an original dynamics 

dXt = {nXt + 47VE(Xt Xt)Xt) dt + V2dBt (19) 

visiting the whole space M , we may consider an approximation of ([5]) as a projection of the 
same original dynamics on a "manifold" defined by the constraint "E||Xt|p constant". The 
difficulty is that giving a mathematical meaning to the latter constraint is not straight- 
forward. The aim of this section is to give a proper meaning to this projection, and to 
identify the projected dynamics with the dynamics (I13p - (jl4p that leads to the Doi closure. 
To keep our exposition simple, we omit the nonlinear term in the drift term of (J19p (namely 
A^ = 0). The reasoning below generalizes to the full drift. 
The approach we propose is to consider / > 1 replicas 

dXi = KXidt + V2dBi 

(for 1 < i < /) of the dynamics (I19p and project the system obtained on the manifold 



tEii^^ 



^"2^L2. 



j=l 



We thus impose that the empirical average is constant, and we are interested in the limit 
/ —7- oo. Of course, the d-dimensional Brownian motions Bl are assumed to be independent. 
The projection is performed using the D'Alembert Principle. Indeed, the constraining 
force does not bring or subtract energy from the system: it is directed orthogonally to 
the submanifold on which the constrained system evolves. More precisely, denoting by 
Xi = {Xl, . . . , Xj) £ M.'^^ , the projected dynamics writes (see for example [U I13|): 

dXt = P(Xi)KXi dt + V2P(Xt)dBt - {dl - l)-^dt, (20) 

where K is the dl x dl block diagonal matrix composed of the blocks k of size d x d, and 
P(X) is still defined by ([3]), with X e R'^^ . We fix 



4=1 



2 _ r2 



L', (21) 



at initial time and this quantity is by construction preserved in time. We also assume that 
the random variables Xq are identically distributed so that, from (j2ip . 



E{\\X^r) = L\ (22) 

As mentioned above. Equation ([5]) is recovered using this projection procedure with only 
one replica: 1 = 1. Here, we consider the limit / — t- oo. 

We now pick the first component X/ € M of our vector Xj and consider its evolution 
equation 

dXl = ^X\ dt + V2dBl - X/^^^-i^ dt - Xl'^^ - {dl - 1)^ 



XII2 t IIY 112 \ ' IIY 112 



Since ||X+|P = I L^, this also writes 



dXl = nXl dt + V2dBl -^\J2XI. kXI dt - ^±^^i^ - ( d - j) ^dt. 

L^ I ^ L y/i ||Xi|| V I) !?■ 

In the Hmit / — )■ oo, we formahy obtain that X\ converges to Yj solution to 

dYt = KYt dt + V2dBl - ^ ¥.{Yt ■ nYt) dt - d ^dt. (23) 

This limit may be rigorously justified as follows. 

Proposition 3.1 Let X[ = {X^ ' , . . . ,X^' ) G M."^^ be a solution to PD|) (we here explicitly 
indicate in the superscript the dependence on the number of replicas I) and Yf solution 
to (j23p . The initial condition Yq is assumed to satisfy 



E{\\Yof) = L^ 
and EdlYoll*) < oo. Consider i.i.d. copies Ig o/Iq and define the initial condition 

I \ -1/2 

so that (Xq )i<i</ are identically distributed random variables satisfying (j2ip . Then, for 
any positive time T > 0, there exists C > such that, for all positive / G N, 





The proof is provided in the appendix. Of course, the convergence result holds for any 
component XI of the vector Xf, and by standard results in propagation of chaos |17| . we 
actually have that any subset of components {Xl^ , . . . X^'') converges (in the limit I — )• oo) 
to (1^*^, . . . ,Y^''), where the processes (1^*) are independent copies of Yt solution to (|23|) . 

Notice that since E(||yt|p) = L^, we may therefore equally well write the dynamics 
on Yt in the following form 

dYt = KYt dt + V2dBl - ^^^IE(y, . nYt) dt - d^^^^^dt. 

This agrees with (J13p for {R = Id, A = d), thereby providing a justification of our particular 
choice p^ in the previous section. 

To summarize, the original model for rigid rods ([5]) may be seen as a projection of the 
dynamics ()19p onto the submanifold ||Xj|p = ||Xo|p, while the approximated model (fT3|) 
for [R = Id, A = d) which is consistent with the Doi closure ([8]) may be seen as the original 
dynamics (J19p constrained to have a fixed average length: E (||Xtp) = E (||Xo|p). This 
yiels a microscopic interpretation of the Doi closure. 

4 Long-time behaviour of our approximate model 

We henceforth consider the model (|13p - (|14p (namely [R = Id, A = d)) which we have 
built from the original model ([5]) by approximation. Throughout this section, we work in 
dimension d = 2. This is a crucial assumption, specifically needed for our technique of 

10 



proof which makes use of the Poincare-Bendixson Theorem. Additionahy, we assume that 
the matrix k is defined by (jl]): 



Pe 
T 



a + 1 
a-l 



and that the initial condition Xq satisfies E(||Xo|p) = L^ = 1 (so that, for all positive 
time, E(||Xi|p) = 1). Given these assumptions, we now recall, for the convenience of our 
reader and the consistency of the present section, the model under study: 

dXt = {nXt - K : E(Xj ® Xt) Xt) dt 

+ m {E{Xt ® Xt)Xt - E{Xt ^ Xt) : E{Xt Xt) Xt) dt (24) 

+ V2dBt-2Xtdt. 

We also recall that the conformation tensor M(t) = E{Xt ^Xt) then satisfies the ordinary 
differential equation: 

-p = {kM + mJ -2k: mm) 

+ 8N{M'^ -M : mm) (^^) 

+ 4(Id/2-M), 

where tr(M(t)) = tr(M(0)) = 1. The (non-linear) Fokker-Planck formulation associated 
to (1241) and established in (1151) writes: 



dtp 

-— = div {{—KX + K : M[ip]x) tp) 

+ 4iVdiv {{-M[7p]x + M[7p] : M[ip]x) ip) (^^) 

+ AV- + 2 div (xV') , 



where 

M\ip{t,-)]= I x®xip{t,x)dx 



and tr(M[V'(0,-)]) = 1- Notice that t ^ M['0(t, •)] then satisfies ([25]). The aim of this 
section is to study the longtime behaviour of the solution ^ to the Fokker-Planck equa- 
tion ([26]) . 

As it is standard for such analysis, we study the longtime behaviour of a solution to 
the Fokker-Planck equation ([26|) . assumed sufficiently regular for our manipulations to be 
valid. We refer for example to [2j for an appropriate functional setting to justify such 
calculations. 

4.1 Long-time convergence of the solution to ( !25l) to a periodic solution 



We first consider the closed ordinary differential equation (I25p on M, momentarily leav- 
ing (|26p aside. We will show that, under some assumptions on the parameters and the 
initial condition M(0), M converges, in the long time limit, to a periodic solution. This is 
an extension of the result |14^ Theorem 5.1]. We also refer to that contribution for a more 
thorough study of the longtime behaviour of the dynamical system, for other regimes of 
the parameters. 

Proposition 4.1 Assume that Pe is sufficiently small, |a| < 1 and N > j^^- Then there 
exists an open subset 0, of the ensemble of positive definite matrices with trace one (Q will 
be made precise in the course of the proof below), such that: 

11 



• There exists a unique periodic-in-time function Mp^^{t), valued in Vt, solution to (|25|) ; 

• For any initial condition M(0) G Vt, the solution M{t) to (I25p converges to Mper{t) 
exponentially fast in the long time, that is: there exist C, A > 0, such that, for all 
t> 0, 

\\M{t) - Mper(t)|| < Cexp(-At). (27) 

The remainder of this section is devoted to the proof of Proposition 14.11 

Proof. 
Step 1: Existence of a periodic-in-time solution 

Using Proposition 12.11 we know that the solution M(t) to ([7]) is symmetric and satisfies 
tr(M(i)) = 1 since this holds true at initial time. Introducing as in Section 12.31 the 
traceless part Q = M — Id/2 of M, we may always write Q, in the two dimensional setting 
we consider, under the form 



Q{t) 



x{t) y{t) 
y{t) -x{t) 



Now the evolution equation (J12p equivalently reads 



dx 



— = -4x (1 - iV + AN(x^ + y^)) + Pe y(l - 2ax), 



dt 

-1 = -4y (1 - iV + 4N{x^ + y^)) + Pe (-x + 

Introducing the polar coordinates 



2ay' 



(28) 



X = rcosip, 
y = r sin 99, 



(29) 



we rewrite (1281) in the form 



-Pe I 1 - — cos (/? ) , 



aPe, 



-4r (A^(4r2 - 1) + l) + (1 - Ar"^) gin ip. 



(30) 



We now consider two positive constants ei and €2 satisfying < ei < -^j^ and < £2 < 



N-l 



jj^ respectively. Set ri = y ^^j^ — ei and r2 = \/ ^-^j^ + £2- Note that, by construction, 
< ri < r2 < 2' Then, if ri < r(0) < r2, one has 

ri < r(t) < r2 

for all positive times, provided that Pe is sufficiently small. Indeed, it is easy to check that 
if Pe < Pe, where 



— IGiVei 

Pe = — r-, — iiiin 



ri 



r2 



i, + 4ei'i,-4e2 



>0, 



N 



N 



then for all angles 93, 



aPe aPe / 1 

-4ri {N{4rf - 1) + l) H (1 -4r^)sin99= 16ri7Vei H ( h 4ei ) s'mip > 



12 



and likewise, 



aPe aPe /I \ 

-4r2 (7V(4r^ - 1) + l) + — -(1 - 4r|) siiK/? = -16r2iVe2 + — - ( — -4e2)sinv9 < 0. 



This shows that — 
dt 



, dr 
> and — 



< 0, thus the annulus 

r=r2 



n 



Ux,y),ri < \J x^ + 2/2 < r2\ 



is stable under the flow (for positive time). The domain $7 mentioned above in the Propo- 
sition [4T] is now made precise and defined by: 



(x,y) G ^^ M 



X y 

y -X 






We now want to show that there is no stationary point in the annulus il. Since, by 
assumption, A^ > . _ ^ , we may assume that ei > is chosen sufficiently small so that 

1 

N > 



1 - 4ei - a? 

which is equivalent to — < 1 in il. This implies that 

2r 

^ = -Pe(l-|:Cosv.)<0 (31) 

since the solution remains in Q, and thus that there is no stationary point in Q. 

From the Poincare-Bendixson Theorem (see for example |15l Theorem 6.12]), we then 
obtain that, for any trajectory with initial condition (3:0,2/0) in ^j its w-limit set is a 
periodic orbit, that is, the trajectory of a periodic solution. We recall, for consistency, that 
a point {xoD,yoo) is in the w-limit set of the (forward) trajectory starting from (xo,yo) if 
there is an increasing sequence of times t„ going to infinity and such that (x{tn),y{tn)) 
converges to (xqo , ?/oo ) as n goes to infinity. 

A corollary of the previous statement is that there exists a periodic solution {xper,yper){t) 
to the equation (128p in il, and thus an associated periodic solution Mper(i) to (|25p . 

Step 2: Properties of the periodic-in-time solution 

In order to prove the long time convergence to the periodic solution and make precise 
the rate of that convergence, we now compute the divergence of the vector field in the 
right-hand-side of (f28]l : for (x,y) G fi, 

Q 

D{x, y) = jr- (-4x (1 - A^ + 4A^(x2 + y"^)) + Pe y(l - 2ax)) 

+ — (-4y [1-N + m{x^ + y^)) + Pe (-x + ^ - 2ay'^ 

= 8 (Af - 1) - UN{x'^ + y2) _ 6Pe ay 
<S{N -I)- MNrl + 6Pe a rs 
= -?>{N -l) + UNei + 6Pe ar2. 

We thus see that if ei and Pe are chosen sufficiently small, then 

D{x,y)<Q 
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in Q. Following |12| I15| . we deduce that the equation (|28p has a unique stable periodic 
orbit in Q. The uniqueness follows from a generalization of the Dulac criterion. Let us 
recall the main arguments for the stability statement. We introduce the Poincare map 
associated to the first return to the section 

S = {{x,y) £ fl such that x > 0,y = 0} 

of Q. It is a standard result |12^ Equation (1.17)] or |15t Equation (4.51)] that 

/) = exp f / D{Xper{t),yper{t)) dt j , (32) 



where T is the period of the periodic solution {xper,yper){t), gives the derivative of the 
Poincare map at its stationary point (namely the point where the periodic orbit intersects 
S). Thus, D < implies that the derivative is strictly smaller than one {p < 1), which 
yields the exponential convergence to the stationary point of the Poincare map, and thus 
the (exponential) asymptotic stability of the periodic orbit. 

For further use (see the proof of Proposition 14. 5p ). we establish uniform bounds on 
the eigenvalues of Mp^r- We first note that Mp^r is a symmetric positive matrix. It is 
a simple consequence of the fact that {xper{t),yperit)) G f^- Indeed, one can check that 
tr(Mper(i)) = 1 and 

det(Mpe.(t)) = ^ - ixl,,{t) + yl,,{t)) > 1 - ri > 0. 



This implies that the eigenvalues of Mper{t) are ^ i \ Xpf,j.{t) + yp^j.{t), and thus bounded 
both from below and from above. In the sense of symmetric matrices, 

< Q - rz") Id < Mp,r{t) <(l + ^2) Id. (33) 

Step 3: Convergence in the long time 

To conclude our proof, we now show the convergence (|27p . Consider a solution {x,y){t) 

to (j28p and the sequence 

{x,y){tk) = (xfc,0) 

of its successive return points to the section S. Otherwise stated, x^+i is the image of x/^ by 
the Poincare map. Notice that x^ is a monotonic sequence (since two trajectories cannot 
cross). As explained in Step 2, we also know that the sequence Xk converges exponentially 
fast to the fixed point x* of the Poincare map: there exists C > and p G (0, 1) (which 
can be chosen arbitrarily close to p) such that, for all k >0, 

\xk-x*\<Cp''. (34) 

Since (x*,0) is on the periodic orbit, there exists a time t* G [0,T) such that 

yXperiUper/yt ) — \X ,Uj. 



Without loss of generality, we may assume t* = 0. We now remark, using (|3ip . that 
if : [tfc,tfc+i) I— 7- [0, 27r) is a one-to-one (actually strictly decreasing) function. We may 
therefore use ip itself to reparameterize the trajectory between two successive return points. 
It follows from dSOl) that 



1 f^^ 1 



tk+l -tk- ^ I -. a n^^- 

Pe Jo l-2K^^T(e))Cose 



14 



Likewise, along the periodic trajectory (rpen ¥'per)(i) (namely {xper,yper){t) in polar coor- 
dinate), we have: 



1 /"ZTr 



^ ~ Pe io 1 r-^7STT cos/^' 

We now use the fact that {r,(f){tk) = {xk,0) is close to {rper,fper){0) = {x*,0) (by virtue 
of (j34p ) and the Lipschitz property of the flow associated to (|30p as a function of the initial 
conditions over a finite time interval, to get that 

\tk+i -tk-T\< Cp'' 

where here and below C > denotes irrelevant constants. 

This implies that there exists a time Tq such that for all A; > 0, 

\tk-To-kT\<Cp''. 

Indeed, ife - to = Eti(*i -ti-i) = kT + Eti(*« -ti-i-T) = kT + E/>i(i« - i«~i - 
^) - Y.i>ki'ti -ti~i-T). Denoting Tq = to + Ez>i(*/ - i«-i - 7"), we have tk-To-kT = 
— '^i^kiti — ti-i — T), thus our claim. 

Then, we have: for k > 0, 

\\{Xper,yper){kT) - (x,y){To + kT)\\ 

< \\ixper,ypeT)ikT) - (x,y)(tfc)|| + \\{x,y){tk) - {x,y){To + kT)\\ 

<\x*-Xk\+C\tk-{n + kT)\ 

<Cp''. 

Using the Lipschitz property of the flow, we thus get: for all /c > 1 and for all t G 

[kT, (k + 1)T) 

\\{Xper,yper){t) - (x, y)(To + t) || < C p'' . 

This implies (I27p with A = — ln(p). 

4.2 Analysis of the Fokker- Planck equation fl26|) for M{t) given 
In this section, we consider (|26p for a given M{t): 

dtp 



div {{—Kx + K, : Mx) ip) 

+ 4Afdiv {{-Mx + M : Mx) tp) 
+ A'ip + 2 div {xip) . 



(35) 



This equation can be rewritten in the form 

^ = dw{K{t)xi;) + AiP (36) 

ot 

with 

K = -K + K : Mid + AN {-M + M : Mid) + 2ld. (37) 

We are thus considering here a linear Fokker-Planck equation. 
We begin with a basic remark, see for example [H [3]. 
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Proposition 4.2 Let ipi and%l)2 he two solutions to the Fokker-Planck equation (|36p . Then 

at 
where 

i:f(V'ilV'2) = y"in(V'i/V'2)V'i 

denotes the relative entropy (of ipi with respect to 1^2), o-fid 

I{^i\i'2)= /|Vln(Vi/V2)|'V'i 



is the Fisher information (of ^i with respect to ■ip2)- 

Proof. Set b{t, x) = K{t)x. We argue formally. Our manipulations are standard and can 
be made rigorous using appropriate functional spaces and cut-off functions. We write 



d 
dt 



H{Mi^2) = j dt^l^i " /^*^2|^ + /^'' (1^) ^'"^^ 

= - /" div(6'02 + VV'2)^ + I div(6^i + V^i) In (^ 
= j{hi^2 + VV'2) • v|i - j{b^ 



H+VVi)-Vln(^ 

^^2 

VV^2-V^- /vVi-Vln( ^ 

W2 J \V2 

^VV^2-VVi) -Vlnf^ 

V2 J \W2 

= -i{ipi\i>2). 

Notice that this proof does not require M to satisfy (|25p . ^ 

We now build an explicit Gaussian solution to (I35p . 



Proposition 4.3 Let M{t) he a given time- dependent symmetric definite positive matrix 

with tr(M(0)) = 1. Introduce the associated two-dimensional Gaussian prohahility density 

function 

, , , v^detM-i(t) / x^M-i(t)x 

V'Af (i, x) = ^ ^ exp ' ^ ' 



27r 'V 2 

Then, ipM satisfies (I35p (for the given function M(t)) if and only if M(t) satisfies (i25|) . 

Proof. Let us denote P{t) = M'^it), so that 



, , , Jdet P(t) f x^P(t)x 
Ztt 



It is easy to calculate that 



x^dtPx , 9t [Vd^tP] / ^Tp^ 



dti'M = 7, V'M H 7, exp 

Z ZTT 
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where here and for the rest of this proof we use the short-hand notation dt instead of 4^. 
Note that (using the symmetry of P) , 

dt (Vdetp] = ]- U/detP\ ^ ^((det P) 

= ^\/detP5i(lndetP) = ^Vdet Ptr(P~^atP). 



It follows that 

„ , , x^{dtP)x , 1 



+ hiip-^dtP)] i^M- 



2 
Plugging ^jvf in the right hand side of the Fokker-Planck equation (j26p yields 

div(ifx^Af + VV-A/) = tr(i<f - P)V'Af - x^PiK - P)xi;M, 
where K is defined by (|37|) . so that ipM is solution to (I26|) if and only if 

_x^_{dtP)x ^ l^^(p-iQ^p) ^ ^^(^ -P)- x'^P{K - P)x, 
for all X G M^ . This is equivalent to the couple of conditions 

-tr(p-^9tP) =tr(i^-P), 

2 ^ * ^ ^ '' (38) 

dtP = PK^K^P -2P'^, 

where, for the second line, we have equated the symmetric part of the two second order 
coefficients. We immediately remark that the second line of (|38p implies the first line, by 
elementary properties of the trace. 
We now write (1251) under the form 



dtM = -{KM + MK'^) + 2Id, 

where, again, K is defined by (I37p . Thus M satisfies (I25p if and only if P = M~^ satisfies 

dtP = PK + K^P - 2P2 , (39) 

where we use the fact dtP = —PdtMP. 

By comparing (p8j) and ([39]) . we thus get the result: Vm satisfies (p6|) if and only if 
M{t) satisfies ([25]). 

Remark 2 Notice that this proposition implies the existence of Gaussian solutions to the 
non-linear Fokker-Planck equation (j26p . Indeed, let us consider a Gaussian initial condition 

m-) and the assoc^ated M(0) = /x^x^(0,.)dx ^^al condMon to J25j. Let M^t) 

he the solution to (j25p with this initial condition. Let us then consider ipM the Gaussian 
solution to ()35p built in the previous proposition. We notice that tpM is then a solution to 
the non-linear Fokker-Planck equation ()26|) . by uniqueness of the solution to ()25p . 



We now proceed with a uniqueness result for the periodic solution to (j26p . This result, 
for which we provide here a self contained proof, is also a consequence of the convergence 
stated in Proposition 14.61 and proved in the next section. 



17 



Proposition 4.4 Consider ipper{t,x) a (sufficiently regular) periodic solution to the non- 
linear Fokker-Planck equation (j26p . Define the associated periodic time-dependent matrix 



Mper(i)=M[V'per(t,-)] 



X (g) Xlljperit-iX) dx. 



Assume that Mper(O) G ^, where Q. has been introduced in Proposition \4-l\ Then, Mp^.^. is 
the unique periodic solution with value in Vt built in Proposition \4.1\ and 



t 



per 



i^KU 



is the associated Gaussian solution built in Proposition \^. 



Proof. The uniqueness of the periodic orbit in Q ensures that Mp^r is indeed the periodic 
solution built in Proposition 14.11 (up to a phase change we may ignore). 

Then, by uniqueness of solutions to (l25l) . we notice that ij^per and ipMper ^^^ both 
solutions to the same linear Fokker-Planck equation (j35p . with M(t) = Mper{t). The 
question is thus now the following: given the function Mp^r, we have to show that the 
periodic solution to (|35p is unique (up to a normalization factor of course) and equal 

to VMper ■ 

Let us denote ipper a periodic solution. We want to show that ipper = ipMper- Arguing 
as in Proposition 14.21 we have 



d_ 
di 



H{ijper\i^Mp,r) + Hi'perllpMp, 



0. 



(40) 



If ^lJper and ipMp^r share the same period, say T, then the function H^ipperli^Mper) i^ also 
T-periodic in time, and thus integrating the above equation from to T, we obtain 



I{'lpper\lpMp,r)dt = 0. 



Since the function / is nonnegative, this immediately implies that tpper = ipMper o^^ [^i^] 
and thus for any time. 

If the period T of ipper is different from the period T of V'Mperi '^'^ slightly adapt the 
above argument. From standard results on continuous fractions (see for example \9\ il6j). 
there exists sequences of integers p„, Qn such that 



In rp Pn 



1 2 

< — < 



ln+1 



n+1 



where 4> 
and linir, 



2 ■ Thus, if we set t„ = PnT^ we have r^ = q-nT + En where lim„^oo £n = 
3o T„ = CO. Then, we have, for n sufficiently large, 

l^(V'perlV'Mper) K) - H{%l)per\^Mp,r) (0)| 

< H{ll)per\tpMp,r) (Tn) - H {^l^per {PnT)\i^ Mp,r (QuT)) 

+ H{i;per{PnT)\4^Mp,r{(lnf)) - H{tpper\i^Mper) (0) 

H{'4)per{PnT)\'4)Mp,MnT + En)) - H{4)per{PnT)\4^Mper{QnT)) 



1pper{0,-)'^T^ 



i^Mperi^n,-) 
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and the right-hand side converges to zero as n goes to infinity, using the continuity in time 
of Mper, and the fact that J'i/'per-(0, x)||x|p dx < c«. 
Integrating (|40p in time from to t„, we thus get 



hm / I{lpper\lljMj,er) dt = 0, 

and thus 

/■max(T,T) 

Jo 

since / is nonnegative and r„ > max(r, T) for large enough n. This shows that ^per = 
V'M e,. on the time interval (0,max(r, T)), and therefore for all times. {> 

We conclude this section with an inequality that will be useful below to show exponen- 
tial convergence to periodic solutions for (I26p . 



Proposition 4.5 Let Mp^r he the periodic solution to (j25p . and Tpp^r = ipMper ^^^ asso- 
ciated (Gaussian) periodic solution to (I26p . Then, ipper satisfies a uniform in time log- 
arithmic Sobolev inequality in the following sense: there exists /u > such that for any 
probability density function ip and for any time t E [0,T), 

H{i;\lljperit, •)) < ^I{^P\^|Jper{t, •))• (41) 

Proof. It is well known that centered Gaussian distributions with covariance matrix M 
satisfy a logarithmic Sobolev inequality with parameter the inverse of the largest eigenvalue 
of M, see for example [1]. 

Now, (|33p precisely shows that the eigenvalues of M are uniformly bounded from above 
by a time-independent constant. This concludes the proof. ^ 

4.3 Long time convergence of the solution to the non-Hnear Fokker- 
Planck equation ( !26|) to a periodic solution 

Our final step is to address the long time behaviour of the solution to the non-linear Fokker- 
Planck equation (I26p . Assume the parameters (a, A^, Pe) are such that the conclusion of 
Proposition 14.11 holds: M{t) = Mliplt, •)] converges exponentially fast to the periodic 
solution Mper{t) (with value in $7) of ([7]). Let us denote 



i^per = V'Mi 



per 



the associated unique periodic solution to (I26p (see Proposition 14.41 above) . Then, consider 
a solution %[) of (|26p and assume that the initial condition V'(O) satisfies 



x(g)x^(0,a;)(ix = M(0) G f^, (42) 

where Vt is defined above. We have: 

Proposition 4.6 Under the assumptions of Proposition [777] (and in particular (I42p ). the 
solution ijj to (j26p (which we assume sufficiently regular) converges exponentially fast to 
tpper in the following sense: there exist C > and v > Q such that, for all time t > 0, 

H{il;{t,-)\'4^per{t,-)) < Cexp{-ut). 
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Proof. To the function tp is associated M(t) = M[ip{t, •)] satisfying ([7]), which (by Propo- 
sition |3?T|) converges exponentially fast to Mp^r'- Vt > 

\\M{t) - Mper(t)|| < Cexp(-At). (43) 

The function ip satisfies the linear Fokker-Planck equation: 

-^ = div(K(i)xVj) + AV', (44) 

where K is defined by (|37p . Likewise, the function ij^per satisfies the linear Fokker-Planck 
equation: 

dijj 

— = div{Kper{t)xil)) + Alp, 

where Kp^r is the periodic function defined by (I37p with M = Mper- Notice that from (I43p . 
we get 

\\K{t) - Kperim < Cexp(-At). (45) 

Now, adapting the proof of Proposition 14.21 and rewriting (144 p as: 

-^ = diviKperit)xij) + AV + div([i^(t) - Kperit)]x^), 

we have, for < £ < 1, 

^i/(^|V'per) = -imi^per) + / div([i^(t) - Kp,r{t)]x^P) In (-^ 

= -/(ViV'per) - [[Kit) - Kperit)]x • V In (-^^ V 

< -(1 - e)/(^|Vper) + ^ / \\[K{t) - Kperit)]xf^P 

< -(1 - e)Iii^\^Pper) + j-jK{t) - Kp,r{t)f J ||x|| V- 

Now, using (j43l) . the fact that f \\x\\'^ip = tr(M) = 1, and Proposition 14.51 we get 

j^H{'iP\iPper) < -^^^y^H{i;\ijper) + ^exp(-2At), 

from which we deduce the exponential convergence of H (iplipp^.r) to zero, using the Gronwall 
Lemma. 

Remark 3 It is easy, by making precise all the constants used in the bounds above, to give 
an expression for the rate of convergence v in terms of p defined in ([32]) and /u defined 
in (1411). 



'per 



A Proof of Proposition 13.11 

We provide in this appendix a proof of Proposition 13.11 Adapting a standard coupling ap- 
proach, see for example [10], we introduce N independent copies of the nonlinear stochastic 
differential equation (|23p : 

dYi = kY; dt + V2dBi - ^ E(F/ • kF/) dt - d ^dt, 
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^ij\ 



driven by the same Brownian motions as the processes {X^^)fyQ. Let us recall that, 
from ()2ip . we have: for all positive time t, 



\Y.\\^i 



iJu2^^2 



(46) 



so that 



i=l 



E||X*'^f = L^ 



since the law of the stochastic process X^ is invariant under permutation of the indices 
of its components (X^' ,... ,X^' ). Moreover, since EHYqIP = L^, we also have, for all 
positive time t, 

E(||y/||2)=E(||y,||2) = L2. (47) 

This originates from the fact that E(||yi|p) solves the ordinary differential equation: 

Imnf) = 2E,r, ^ .r,) (i - mp) + 2. (i - « 



We now introduce the difference 



Zl:' = X];' - 1-, 



with initial condition Zq =Yq IL ( j X]j=i \\Y^ 



i\\2 



-1/2 



1 I . We have: 



yi,I 



^0 +[^-Jy] I ^sds 



12 /* ( ^'s'j E ^s' ■ ^xi'' - y: e{y: . ky:) j ds 



,0 L V/ IIX^II IL^Jo 

Thus, using the Doob inequality to get 



xyds. 



E ( sup / ||X 

.s<t JO 



1,1 U ^r -dBr 



ix^l 



<CE( nixyi^^'-'"^^' 



ll^rl 



C I E\\Xyfdr<CL'^T, 
'o 



(where here and throughout this proof, C denotes irrelevant constants that depend on T, 
K, d, L but not on /), we obtain, for any time t G [0,T], (T is fixed). 



,-f||2 



s<t 



i,I\\2 



E sup||Z*''f <E||Zo''f + C7 / E sup||Z, 



:,-fii2 



t 

+ C I E 
/o 



\r< 

I 



ds 



^F J E ^s' ■ ^^s' - ^s ny: ■ ^i: 



i=i 



ds (48) 



+ 






We are now going to estimate from above the first and third terms of the right-hand side. 
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^i.I 



We begin with the first term, which involves the initial condition Zq . Using that, for 
all x,y > 0, (x — 2/)^ < \x'^ — y^|, we have: 



milXg II — ill I 11^0 I 



L 



1 Y^-f ||V«II2 
/ 1^1=1 IMO II 



E yEll^ol 



L 



i=i 



E 



E 



1 Y^J \\vi\ 
I 1^1=1 iMol 

2N 



(^-i^)) 



r2 1 v^/ ||v« 



L+j'jEUwy^ 



^i\\2 



< 



^AW--,T.p-if 



i=l 



lVar(p^ 
/ L2 



We now consider the term (l48ll . which we split as follows: 
1 ^ 1 ^ 

x'/j Y. xi'^ ■ KXi'' - y: e(y: . ky:) = (xf - y:)- ^ xi^' • kx^'' 



i=i 



i=i 



1 ^ 



i=i 



+ ^^M 7 E ^^ '^^/' - K(^^ ^^^') ) 



For the first term (|50p . we have, using (jl6 



E 



1 I 
{Xi:' -Yi)-Y,Xj:' ■ kX^' 



i=i 



:,-fii2 



< CEllZ^f < CE sup||Z, 



'J||2 



r<s 



For the third term (1521). we write: 



E 



(49) 



(50) 
(51) 
(52) 



(53) 



= 72 E ^^ (r;f (y/ • ^y/ - E(y/ • ky^)) {y^ • ky^ - ^Ys ■ ^y' 
j,k=l 

By independence of the stochastic processes (l^*)j>i, the terms in the sum are zero if j ^ k. 
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Thus, 



E 



y: f jiz^s' ■ ^Yi - E(y; • ky:) ) = -1 ^E (||y;f (y/ • ^y/ - E(y/ • nYi)f) 



< 72 E V^(r;r)viE (i^/ • '^n-'' - E(y/ • KYi\ 



< 



7 



(54) 



using that supg^rg j^i E (||ys||^) < oo, which is easy to check from (j23p provided the initial 
condition Yq is assumed to have finite moments up to order 8. We finally estimate the 
second term (1511). We have 



E Y:\Y.^s'-^x^'-yi-^i] 

< CE I ||y;f 1 r£\\zi''\\ {wxi^'w + ||y/||) J J 



< CE 



C 
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HEJ=ir/lP<2L2 



||Y'i||2>2L2 



Consider the second term (the first term is addressed similarly). We write 



|2 I \\VJ\\'^\ 1 1 

+ 11-^. II ) I ^iE-=iri'IP>2i2 






^il|y/||2>2L2 

3 

=il|y/lP>2L2 



<CE(sup||Z^'^||2')+^. (55) 



Here, we have used a concentration inequality, based on (l47|l : for n = 2, 3 



< 2""^E 



+ 2" ^L^ng^^^^^^^^^^^^^^^^^^^^^ 



C 



The last line is a consequence of the following: for ^j = WYsW"^ — L^ i.i.d. centered random 
variables with finite fourth moment, a = L^ > 0, and m = 0, 2, 3, 



'{{ip) HEt..>«)^^-(()p.)) 



4\ 

<72' 



where C depends on sup^gro^rj E ( ||y/||'^ j < oo. 
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Inserting l^-l^-l^-l^ in (gS]), we obtain: for all time t G [0, T] 
E(sup||ZJ'^in <C [ E(sup||Z;'^|n ds + ^, 



\s<t / Jo \r<s J ' 

and we conclude using the Gronwall lemma. 
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